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ABSTRACT 

A flat Universe model supported by recent observations has 18 possible choices for its overall topology. To detect or exclude these 
possibilities is one of the most important tasks in modern cosmology, but it has been very difficult for globally inhomogeneous 
ones because of a long-time calculation. In this brief paper we provide an object-based 3D method to overcome the problem, as an 
extension of Fujii & Yoshii (201 la). Though the test depends on the observer's location in the universe, this method drastically reduces 
calculation times to constrain inhomogeneous topologies, and will be useful in exhaustively constraining the size of the Universe. 
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1. Introduction 

Theory of modern cosmology is based on Einstein's General 
Relativity, and recent observations favor a A-CDM cosmol- 
ogy with vanishing curvature (e.g., Q tot = 1 . 0050^^0 from 
WMAP+BAO+SN data, by Hinshaw et al. 2009), which suc- 
cessfully describes the observed properties such as the cosmic 
structure formation, the cosmic microwave background (CMB) 
anisotropies, and the accelerating expansion. However, General 
Relativity does not distinguish two spaces with the same curva- 
ture but with different topologies. Though the curvature of our 
Universe is exactly zero, we still have 18 choices for the overall 
topology of the Universe (Nowacki 1934), e.g., the multicon- 
nected 3 -torus T 3 with finite volume. 

A multiconnected space is a quotient space of the simply 
connected space with the same curvature (§ 3 ,E 3 , or H 3 ), by a 
holonomy group T. An observer in the space would see it as a 
finite or infinite 2/T-polyhedron (the Dirichlet domain) whose 
K pairs of faces are glued by the holonomies. There exist the 
second, third, and more, geodesies between a given object and 
the observer, which across the Dirichlet domain. As a result, 
multiple images of single objects, often referred to as "ghosts", 
can be observed in a multiconnected space (for details, see, e.g., 
Lachieze-Rey and Luminet 1995). 

Based on this prediction, many object-based works for con- 
straining cosmic topology were carried out, i.e., they searched 
for periodic and symmetrical patterns made by ghosts of galax- 
ies, galaxy clusters, or active galactic nuclei (e.g., Fagundes & 
Wichoski 1987; Demiariski & Lapucha 1987; Lehoucq et al. 
1996; Roukema 1996; Uzan et al. 1999; Weatherley et al. 2003; 
Marecki et al. 2005; Menzies & Mathews 2005). However, their 
methods are valid only for particular topologies and/or, not sen- 
sitive enough to constrain spaces that are comparable to the ob- 
served region in size. The method by Roukema (1996) has both 
the generality and the high sensitivity, but the algorithm is not so 
sophisticated and the calculations take very long times. 



One of the recent trends is to use the circles-in-the-sky 
method (Cornish et al. 1998) that is to search for intersections 
of the last-scattering surface and the faces of our Dirichlet do- 
main. They appear as circles with the same temperature fluctua- 
tion pattern in the CMB map since they are physically identical. 
This method can detect any topologies, but checking all the pos- 
sibilities requires an extremely long-time calculation because of 
the many free parameters: radius of the matched circles and ce- 
lestial positions of centers of the two circles. Such exhaustive 
analyzations have not carried out yet, and various authors have 
searched for matched circles limiting to antipodal or nearly an- 
tipodal ones, using the WMAP satellite's data. However, they ob- 
tained diverse results (e.g., Aurich 2008; Roukema et al. 2008; 
Cornish et al. 2004; Key et al. 2007; Bielewicz & Banday 201 1), 
which suggests that there are methodological problems. 

One way to break down this situation is to use observational 
data that is independent of the CMB observations, e.g., astro- 
nomical objects such as galaxies and quasars. This motivates 
us to revisit the object-based 3D methods, and we developed a 
method that can constrain any of the 18 flat topologies with sim- 
ple algorithms, and is much more sensitive to topological signa- 
tures than the preceding ones (Fujii & Yoshii 2011a, hereafter 
FY1 la). Nevertheless, it also has a similar problem in checking 
all the possible topologies: calculation time becomes very long. 
The aim of this brief paper is to provide a technique to overcome 
this problem. Throughout the paper we consider flat universes. 



2. Method 

2.1. Mathematical background and definitions 

Mathematical classification of the holonomy groups for flat 
spaces was completed by Nowacki (1934). Any holonomy y can 
be written as y — yjym, where yj is a parallel translation and 
ym is one of an identity, n-th turn rotations (for n = 2,3,4, 
and 6), and a reflection. Those spaces whose holonomy groups 
include only translations (y#r - id) are called (globally) homo- 
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geneous, otherwise called (globally) inhomogeneous. We inves- 
tigate the latter topologies in this paper. 

We borrow a convenient way of writing the holonomies that 
is to use a 4D coordinate system (w,x,y,z) where the simply 
connected 3-Euclidean space E 3 is represented as a hyperplane 
w — 1, so that a usual 3D vector (x,y,z) is represented as a 4D 
vector* = (l,x,y,z). Every holonomyy is also written as a4x4 
matrix, e.g., a half-turn corkscrew motion is written as: 



y = U 
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where U is a 4 x 4 matrix representing the choices of the coordi- 
nate systems, which reduces to U — id if we choose our z-axis to 
be parallel to the rotational axis (hereafter, we call this straight 
line the fundamental axis of the holonomy, denoted by l/ m ). The 
unit vector rtf m that is parallel to lj m is called the fundamental 
vector, we identify two antipodal vectors n/ m and -«/«„. For a 
holonomy of reflection, fundamental axis lf un is parallel to the 
reflectional plane. 

In the rest of this paper, unless otherwise stated we continue 
to consider the same half-turn corkscrew motion as an example. 

2.2. Summary of the 3D method of FY1 1a 

In this section we review the FYlla's object-based 3D method 
for detecting cosmic topology (see the paper for details). Our 
assumption is that the Universe has the spatial section with zero 
curvature, as suggested by recent observations. This assumption 
is also favored by theory of the quantum creation of the universe 
(Linde 2004). 

If a pair of comoving objects x \ and x\ are linked by a holon- 
omy y, we have by definition 

x\ = yxi = jtJntX] = JntX\ + L, (1) 

where L = (\,L\,La,Lji) is the translational vector. Detecting 
such topological twins requires the parameter searching for 5 
parameters: the translational vector L and the fundamental vec- 
tor n To eliminate L we search for two pairs of ghosts (called 
a topological quadruplet) [(x\, X2), (x\ , x' 2 )] such that 

x\ =y*i= JtJntXx = JntX\ + L, (2) 

x' 2 = yx 2 = y T y NT xi = ym*2 + L. (3) 
Such a quadruplet always satisfies 

x' 2 -x\ = y NT (x2 -*t), (4) 

independent of L. If our z-axis is parallel to the fundamental axis 
lfun (hence U = id), the following relations hold: 

x' 2 - x\ = -(x 2 - xi), (5) 

y' 2 -y\ = -Cv2 -yd, (6) 

z' 2 -z\ = zi -zi, (7) 

since y is a half-turn corkscrew motion here. 

The FYlla's scheme is to search for quadruplets 
[(xi,Xj), (Xk,xi)] that simultaneously satisfy the following three 
conditions: 

1. separation condition: \x/ - Xj\ - \xt - x/\. 



This condition is common to all holonomies, as holonomies 
are isometries that preserve distance. Preceding works such 
as Roukema (1996) and Uzan et al. (1999) also used this 
mathematical property of holonomies. 

2. vectorial condition: eq. (5)-(7). 

This is for a half-turn corkscrew motion; for other types of 
holonomies, see FY 11a. Similar condition for a translation 
was used in Marecki et al. (2005). 

3. lifetime condition: |f,- - \tj - t[\ < tm e . 

The variables ?,-, tj, ■ ■ ■ are cosmic times of objects Xj, 
respectively, and fuf e is the typical lifetime of objects. This 
condition is important when considering short-lived objects, 
e.g., active galactic nuclei. 

Due to these multi-filters, only < 10 pairs of ghosts among the 
total > 1000 objects are sufficient for this test. However, eq. <(5j— 
(0 hold only for the specific coordinate system, so we still have 
to perform the parameter searching for 2 parameters, the funda- 
mental vector «/„„. To our dismay, this requires over about 10 6 
trials that takes a very long time. These strong and weak points 
of the FYlla's method is seen in Fig[TJ where s,- is the num- 
ber of quadruplets that include the object *, as its members and 
satisfy all the conditions 1-3 (see FYlla for details). When us- 
ing the z-axis parallel to «/„„, there can be seen some bumps in 
the histograms that are constituted by topological ghosts, while 
such bumps are not seen when the z-axis deviates from n j u „. The 
property of the simulated catalog used here is described in sec- 
tion 3.1. 




2 4 6 8 10 12 

Si 

Fig. 1. Solid line: the result using the z-axis parallel to «/«„. 
Broken line: the result using the z-axis that is at the angle of 
30° with Hfun. The vertical scale is linear from to 1 and loga- 
rithmic from 1 to 1000. The observer stands at the distance of 3 
Gpc from lj lm (see section 3.2). 



2.3. Fundamental-vector-searching method 

One may want to have a coordinate-independent filter to ex- 
tract topological quadruplets, but it is found to be impossible. 
Consider an arbitrary quadruplet [(jc,-, xj), (Xk,xf)] that already 
satisfies the separation condition and the lifetime condition. It 
then also satisfies the vectorial condition eq. ©-O if we choose 
the z-axis parallel to the vector a + b, where a = Xj — Xj and 
b - Xi - Xu (Fig |2]l. The unit vector n pec oc a + b is called the 
peculiar vector of the quadruplet, for a half-turn corkscrew mo- 
tion. 
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Fig. 2. Even a nontopological quadruplet seems to be linked by 
a half-turn corkscrew motion when we choose the z-axis parallel 
to its peculiar axis n pec oc a + b. 

Therefore it is impossible to judge whether a given quadru- 
plet is topological, however, if there really exists the fundamen- 
tal axis (and vector) in the Universe, the number of quadru- 
plets whose peculiar vectors are parallel to it should be larger 
than stochastically expected. This is the essential idea of the 
coordinate-independent technique described here. Detailed pro- 
cedure, hereafter we call the "fundamental-vector-searching 
(FVS)" method, is as follows. 

1 . The celestial sphere is divided into 12A^ s 2 ide equal area pixels. 
Since we identify two antipodal vectors, the number of pixels 
used is 6N^ ide . 

2. The quadruplets that satisfy both the separation condition 
and the lifetime condition are selected. 

3. For each selected quadruplet, its peculiar vector n pec and the 
celestial pixel containing n pec within it is calculated. 

4. For each object jt,-, the pixels containing more than s m [ n pe- 
culiar axes of the quadruplets that include x, are flagged. 

5. Each pixel is assigned an integer, the number of times of 
being flagged in the previous step. 

These steps practically correspond to counting the number of 
objects with s; > s mm for all the possible z-axes. Note that the 
choices of the parameters N s a e and s m i n are not independent; 
a small N s id e corresponds to large pixels and large stochastic 
noises, so s m \ n should be chosen to be large enough. HEALPix 
scheme (Gorski et al. 2005) is used in pixelization of the celestial 
sphere. 

We consider the half-turn corkscrew motion as an example 
throughout this paper, but the treatments on the other types of 
holonomies are similar except for the step 3. The peculiar axis 
of the quadruplet Xj), (Xk,x{)] for each type of holonomies 
can be calculated as follows. 

1. n-th turn corkscrew motions 

In this case the fundamental axis is defined to be parallel to 
the rotational axis. The apparent angle between a and b, there- 
fore, should be 2n/n when seen from its peculiar vector n pec . 
Such a vector is obtained by rotating the unit vector parallel to 
a + b around the axis parallel to a - b by the angle of 

/ tan0\ 

a — + arcsin , (8) 

I tan II 

where <p is half the intrinsic angle between a and b. The quadru- 
plet with (f> > njn has no peculiar axes, which never occurs when 
n — 2 (half-turn corkscrew motion). 

2. glide reflections 

In this case the fundamental axis is defined to be perpendic- 
ular to the reflectional plane, hence the peculiar vector is the unit 
vector parallel to a - b. 



3. Simulations and Discussions 

3. 1 . Details of simulated catalog 

We generated a full-sky catalog of toy quasars to test the FVS 
method. The properties of the catalog is given here. 

- The standard A-CDM cosmology (Q m = 0.27, Q A = 
0.73, Hq-1\ km/sec/Mpc) is used, so that the local geome- 
try is Euclidean. 

- A half-turn space topology is assumed, such that the 
detectable holonomies are only two half-turn corkscrew 
motions. The translational distance of these detectable 
holonomies is L — 10 Gpc. 

- the comoving number density of quasars obeys the empirical 
law for the luminous quasars (Fan et al. 2001b). 

- Quasar luminosity evolution is simplified such that they emit 
radiation with constant luminosity during the fixed duration 
flife = 10 8 yr. 

- The peculiar motion of quasars is simplified to move with 
constant speed v = 500 km/sec and with randomly chosen 
directions. 

- Other natures of quasars (such as clustering, activity cycle, 
and anisotropic morphology) and the technical uncertainties 
are all ignored. 

We acknowledge that these simple assumptions are not allowed 
for practical applications. Simulations in a more realistic situa- 
tion will be a subject of the coming paper. 

3.2. Main results: detection of the fundamental axis 

First we applied the FY 1 la's method to a toy quasar catalog seen 
from an observer standing at a comoving distance of 3 Gpc from 
lf m . The catalog contains 1010 objects at redshifts z > 4.7, in- 
cluding 12 pairs of ghosts. Two choices of coordinate systems, 
9 = 0° and 30° where 9 is the angle between the chosen z-axis 
and «/«„, were used. These results are given in FigQ]and already 
discussed; topological signal appears only for the former case 
(9 = 0°), and we have to change the z-axis as free parameters 
to detect globally inhomogeneous topologies (see also section 
5.2.1 of FY1 la). This takes an extremely long time, and it is our 
motivation for this work. 

Next we applied the FVS method to the same catalog using 
A^side = 500 and s m j n = 2. If there is a topological quadruplet, 
its peculiar vector n pec ideally coincides with the fundamental 
vector ftf lm , but not in practice due to peculiar motions. The res- 
olution parameter N^de is determined so as to cover this posi- 
tional deviations. The cutoff parameter s m i n = 2 is a reasonable 
choice because the quadruplets with s, > 2 seem to rarely occur 
by chance as can be seen in FigQ] 

The results were smoothed and contoured in 2D maps (Fig 
|3}, where the half celestial sphere is orthogonally projected on 
the xy-plane. As expected, a pixel containing the fundamental 
axis has larger counts than the other pixels whose counts are 
purely stochastic, independently of the coordinate systems. 

This single calculation takes about a few ten minutes with 
a present-day ordinary personal computer, and is equivalent to 
performing about ~ 10 6 trials of the FYlla's method, taking 
about one year. This is the notable advantage of the FVS method 
over the other methods, e.g., the CMB-based circles-in-the-sky 
method that cannot constrain all the possible topologies in a rea- 
sonable time. 
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Fig. 3. The fundamental vector n f u „ appears as a sharp peak in 
both cases; left: 9 = 0°, right: 6 = 30°. The observer stands 
at the distance of 3 Gpc from the fundamental axis lf u „. Plots 
outside the unit disk are fake made by square-shaped meshing of 
the region [-1, 1] x [-1, 1]. 



3.3. Limits of the method: observer's location 

The results of the previous section are based on the somewhat 
idealized assumptions. In practice, there are various effects to be 
considered correctly: physical properties of astronomical objects 
and observational uncertainties. These effects will be treated in 
the coming papers and not here, however, there is another im- 
portant factor that affects this test; the limit of the FVS method 
depends on the observer's location in the universe. 

To see this effect, let us consider a topological quadruplet 
[(*;, X:), (x' : , jc',)], which satisfies 

( J 



a± = -b ± + £_l, 
an = bu - en. 



(9) 
(10) 



where a = x i-X/, b = x'- 



.v,, „ - ~ j .v ; , and s represents the deviation due 
to peculiar motions. The lower indices _L and || represent the vec- 
tor components perpendicular and parallel to lf un , respectively. 
The peculiar vector of this quadruplet is 



oc a + b oc rif un + 



gj| + g-L 

2|fl|,| 



(11) 



If the observer stands near the fundamental axis lf un , the faces 
of his Dirichlet domain are nearly perpendicular to lf m . On the 
other hand, ghosts of the short-lived objects such as quasars and 
starburst galaxies, which are luminous and suitable for high- 
redshift observations, appear near the faces of the Dirichlet do- 
main. As a result, a\\ becomes small and slight changes in po- 
sition due to peculiar velocity induce a large deviation of n pec 
from n fun . 

Therefore, if the Milky Way is situated near //„„, the FVS 
method can no longer detect the topological signatures and a 
troublesome parameter searching is needed. For example, we 
generated a toy quasar catalog seen from the observer standing 
along //„„. The catalog contains 1040 objects including 14 pairs 
of ghosts. In Fig [4] it can be seen that the FVS method does 
not detect the fundamental vector, while the FY 1 la's method 
successfully detects topological signature though in practice a 
long-time calculation to search for n f un is needed. 

Now we discuss the limitation of the FVS method a little 
more in detail. In order that there be a sharp peak, both the fol- 
lowing conditions must be satisfied: 

1. Sufficient number of topological ghosts are in the catalog. 

2. All the peculiar vectors n pec of the topological quadruplets 
are lying in one (or a few) pixel. 
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Fig. 4. In the case where the observer stands along //„„, the FVS 
method does not work well (right), and a tremendous parameter 
searching of the FY1 la's method is necessary (left). 

To satisfy the condition 2, the pixel size should be larger than the 
typical deviation from ii/ B „ of n pec 's of topological quadruplets, 
which can be estimated from eq. (fTTT i. Then the typical number 
of wife's of nontopological quadruplets (including a given object 
jc,) contained in one pixel can be calculated. The condition 1 
claims that this stochastic noise must be sufficiently smaller than 
the number of topological ghost pairs. 

Fig |5] shows the rough estimation of these quantities as a 
function of r b s , the observer's distance from / /„„, in the situation 
described in section 3.1. The typical size of the universe is fixed 
to L — 10 Gpc for simplicity, though it should be treated as a 
variable. The FVS method is valid as long as the stochastic noise 
is small enough compared to the ghost pairs, i.e., 1 Gpc < r b s < 
6 Gpc, which is consistent with the previous results (r b s = 3 Gpc 
and Gpc). When there is no sharp peaks as in Fig |4] we can 
conclude that r Q b s < 1 Gpc or r D b s > 6 Gpc for L — 10 Gpc, and 
the constraints on r \, s for other values of L also can be obtained 
similarly. 
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Fig. 5. The FVS method is not valid if stochastic noise > ghost 
pairs. The universe's typical size is assumed to be L — 10 Gpc 
here. 



3.4. Prospects 

We introduced the FVS method that, along with the FYlla's 
method, is useful for constraining topology of the universe with 
zero curvature. Especially, the former method will provide the 
first observational constraints on the inhomogeneous topologies, 
which take an unrealistically long time when using other CMB- 
based or object-based methods. 

However, for practical application we need more realistic 
simulations. Then our next work will focus on these issues: 
detailed treatments on physical properties of astronomical ob- 
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jects and technical uncertainties of observations. When this ac- 
complished, we will be able to apply the FYlla's method and 
the FVS method to the observed catalogs of our real Universe. 
Recent and future large-scale survey projects, such as 2dF- 
GRS (Croom et al. 2004), SDSS (Schneider et al. 2010), Large 
Synoptic Survey Telescope (LSST) and the Joint Astrophysics 
Nascent Universe Satellite (JANUS), provide sufficient data for 
this test. Formidable progress in techniques for measuring pho- 
tometric redshifts suggests that the spectroscopic surveys are not 
necessarily needed. 
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